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Abstract 

Non-spherical systems described by MOND theories of modified 
gravity arising from generalizations of the Poisson equations are af- 
fected by a MONDian extra-quadrupolar potential 0m even if they are 
isolated and they are in deep Newtonian regime. In general MOND 
theories quickly approaching Newtonian dynamics for accelerations be- 
yond Aqi 4>m is proportional to a multiplicative scaling coefficient a ~ 
1, while in MOND models becoming Newtonian beyond kAq, k ^> 1, 
it is enhanced by k 2 . We analytically work out some orbital effects due 
to <pM in the framework of QUMOND, and compare them with the lat- 
est observational determinations of Solar System's planetary dynamics, 
exoplanets and double lined spectroscopic binary stars. The current 
admissible range for the anomalous perihelion precession of Saturn 
yields \k\ < 2.5 x 10 5 , while the radial velocity of a Cen AB allows to 
infer \k\ < 6.2 x 10 4 (A) and \n\ < 4.2 x 10 4 (B). In evaluating such 
preliminary constraints it must be recalled that QUMOND is not the 
nonrelativistic limit of TeVeS. 



1 Introduction 

The MOdified Newtonian Dynamics (MOND) [1 1 is a theoretical framework 
proposed by Milgrom [2H1] to suitably modify the laws of the gravitational 
interaction in order to explain the observed kinematics of certain astrophys- 
ical extended systems such as various kinds of galaxies [SHZ] and clusters 
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of galaxies [8j. Indeed, their behaviour does not agree with the predictions 
made on the basis of the usual Newtonian inverse-square law applied to 
the electromagnetically detected baryonic matter whose quantity is insuffi- 
cient to accommodate it. In almoslQ all its relativistic formulations, MOND 
implies a sing le acceleration scale [S] A = (1.2 ± 0.27) x 1(T 10 m s" 2 at 
which the laws of gravitation would suffer notable modifications in order 
to explain the aforementioned discrepancy mimicking the effect of the addi- 
tional non-baryonic Dark Matter which is invoked to explain the observed 
discrepancy [10] within the standard theoretical framework. 

Basically, if <?n denotes the Newtonian acceleration within an isolated 
system of total mass M to tai) the modified acceleration g according to MOND 
is 

ffN . 9 m 

where fi(q) is an interpolating function such as 

p{q) — > q for q <C 1 and p{q) — > 1 for q S> 1. (2) 
Actually, the fact that 

— -1^0 (3) 

9N 

H - 1 (4) 



when 
for 



q — > oo (5) 

holds exactly only in some formulations of MOND, such as the "modified 
inertia" ones [llj . and for systems whose matter density p is spherically 
symmetric. Instead, in formulations based on extensions of the Poisson 
equations, such as the nonlinear Poisson model by Bekenstein and Milgrom 
[12] and QUMOND [13] , Milgrom [H] recently showed that g/g^ — 1 remains 
finite in aspherical systems, no matter how ju — 1 becomes small for high q. 
In other words, the dynamics of an isolated, non-spherical system of total 
mass M to t and extension R ^ dyi, wher^l the MOND radius is defined as 



is modified with respect to the Newtonian case even if it is in the deep 
Newtonian regime. 



1 See below for a discussion of this point. 
2 G is the Newtonian constant of gravitation. 
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More precisely, Milgrom |14j considered the case of an aspherical dis- 
tribution of nonrelativistic masses having small extension with respect to 
dyi and whose temporal changes occur on timescales much longer than 
*M = dyi/c, where c is the speed of light in vacuum. By assuming /j, = 1 
to the desired accuracy everywhere within the system, Milgrom [J3], using 
QUMOND [13] , which is the nonrelativistic limit of a certain formulation of 
bimetric MOND (BIMOND) [15J, obtained a qudrupolar correction 4>m to 
the Newtonian potential 

aG ■ ■ 

M (x) = — r iVQij(x) (7) 

with 

Qij(x) = - J p (x'^j (x' 2 Sij - 3x-x^ d 3 x; (8) 

x l ,i = 1,2,3 in eq. ([7]) are the components of the position vector x = 
{x, y, z} of a generic point P in the field with respect to the Center of Mass 
(CM) of the system, while Xj,j = 1,2,3 in eq. ([8]) determine the position 
of the system's mass elements with respect to the CM. The coefficient 
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a= 4 



oo 

w 3/2 \v(w) - 1] dw > (9) 







entering eq. ([7]) depends on a new interpolating function [T 



Hw),w = f (10) 

such as [H] 

v(w) — > w^ 1 / 2 for w — > and v{w) — > 1 for w — > oo; (11) 
i/(w) is related to the MOND interpolating function fi(q) by |14j 

"W-^- (12) 

Milgrom [14] quoted various possible values for a ranging from ~ 0.075 
to pQ ~ 37. Milgrom [TJ] noticed that, in general, it is not necessary to 
have a ~ 1. Indeed, in the original versions of TeVeS [16], Aq is not the 
only characteristic constant appearing in them; in the nonrelativistic limit 
of such theories, the Newtonian regime is obtained in terms of a second, 
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dimensionless constant k 3> 1 when ~ Aq = kAq. Bekenstein and 
Magueijo [17] found k = 1.75 x 10 5 , while Magueijo and Mozaffari [18] used 
k > 1.6 x 10 6 . According to Milgrom |14| . the QUMOND versions of theories 
with k ^> 1 require the use of an interpolating function 



However, the nonrelativistic limit of TeVes-like theories is not QUMOND 
|14j . Finally, we remark that Milgrom [14] felt that theories with n 1 
cannot be considered as generic MOND results. 

As stressed by Milgrom [14], the quadrupolar MOND effect of eq. ([7]) 
considered here has not to be confused with some other MONDian features 
occurring in strong acceleration regime previously examined in literature. 
In particular, it is not the quadrupolar effect [TjJJ[2U] due to the External 
Field Effect (EFE) [2,12, 21] arising when the system under consideration is 
immersed in an external background field; indeed, here the system is consid- 
ered isolated. Even so, residual MONDian effects in the strong acceleration 
regime exist, in general, because of the remaining departure of fi(q) from 

1 when g > 1; their consequences on orbital motions of Solar System ob- 
jects were treated in, e.g., [2l|l9,22,23j. Nonetheless, they are different from 
the presently studied effect, for which it was posed fi = 1 to the desired 
accuracy. Finally, Milgrom [14] showed that the impact of the zero-gravity 
points [17,18,24] existing in high acceleration regions on the dynamics of the 
mass sources themselves is negligible with respect to the effect considered 
here. 

The plan of the paper is as follows. In Section [2] we analytically work 
out some orbital effects caused by eq. ((7|), with eq. (fl4l) . in an isolated two- 
body system. In Section [3] our results are compared to latest observations 
on Solar System planetary motions, extrasolar planets, and spectroscopic 
binary stars. Section U] is devoted to summarizing our findings. 

2 Calculation of some orbital effects 

In the case of a typical non-spherical system such as a localized binary made 
of two point masses M and m with = M tot = M + m, its mass density 

3 It is so because the original versions of TeVeS are not compatible with general rela- 
tivity: they do not reduce to it in the formal limit Aa — > 0. 




(13) 



In this case 



a K = k a K= i, a K= i 



~ 1. 



(14) 
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p can be posed 

p(x) =Mi 3 (x-x M ) + mi 3 (x-x m ), (15) 

where x m and xm denote the position vectors of m and M with respect to 
the system's CM. By calculating eq. ([8]) with eq. (fT5|h eq. ([7]) becomes 

OtG fJj\y r o o 2 2 2 

0M = — 7-5— - 2y m y + 4y M y m y - 

- 2y|/y 2 + x 2 M z 2 + y^z 2 - 2y M y m z 2 + yf f z 2 + (y 2 + z 2 ) - 

- 2x M x m (y 2 + z 2 ) - Qy m yz m z + 6y M yz m z + z„y 2 - 2z^z 2 + 



+ 



(VM ~ ymf + {zm - z m f - 2 (x M ~ x m ) 2 



x 2 - 



- 2 (-3y m yz + 3y M yz + z m y - 2z m z ) z M + (y -2z) z M + 

+ 6 (x m - x M ) [{vm -y m )y + (zm ~ z m ) z] x} . (16) 
where [i\> = mMM^ 1 is the binary's reduced mass. By taking the gradient 



of eq. (1161) with respect to x it is possible to obtain the extra-acceleration 
Ap experienced by an unit mass at a generic point P in the field. The 
extra-accelerations A m and Am felt by m and M can be obtained by cal- 
culating Ap in x m = MM^ 1 r and xm = —mMrr, respectively, where 
f = {x m — XM,ym — yM-,z m — zm} is the relative position vector directed 
from Mtom. We have 



2aA mM 2 

m ~ M b 3< rr ' M 



2aA Mm 2 2 
lM = M 3 ri 3 r ^ (18) 



It can be noticed from eq. ()T7|) -eq. (]T8jl that the extra- forces i^m , i* 1 m acting 
on both m and M are attractive, and their common intensity is 

-- 2aA)/ ^ 3 , (19) 



M h dl 
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in agreement with the result by Milgrom p3]. Finally, the relative extra- 
acceleration is 

A = A m -A M = - 2 -^(j£-)r*r. (20) 



For the following developments, it is useful to remark that, formally, eq. 
(|20p can be derived from the effective potential 

2.1 The pericenter precession 

The longitude of pericenter w = fi+u is a "broken" angle since the longitude 
of the ascending node Q lies in the reference {x, y} plane from the reference 
x direction to the line of the no defl while the argument of pericenter oj 
reckons the position of the point of closest approach in the orbital plane with 
respect to the line of the nodes. The angle w is usually adopted in Solar 
System studies to put constraints on putative modifications of standard 
Newtonian/Einsteinian dynamics [25J. Its Lagrange perturbation equation 

is m\ 



dw\ 1 



dt I n^d 



2 



VT^ \ d{U pert ) tan (I) expert) 
e 8e dl 



(22) 



where U pei t is a small correction to the Newtonian potential; a is the relative 
semimajor axis, e is the orbital eccentricity, and / is the inclination of the 
orbital plane to the reference {x, y} plane. The brackets (. . .) in eq. (|22p de- 
note the average over one full orbital period P h = Ittu^ 1 = 2vry / 'cfiG^M' 1 . 
By adopting eq. (|2T|) as perturbing potential t7 per t in eq. ([22]) . one gets 

As a cross-check of the validity of our result, we repeated the calculation 
of the long-term precession of w by using eq. ()20p as perturbing acceleration 
and the Gauss equations for the variations of the elements: we re-obtained 
eq. ([23]). 



It is the intersection of the orbital plane with the reference {x, y} plane. 
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2.2 The timing in binary radipulsars 



The basic observable in binary pulsar systems is the change r p in the time 
of arrivals (TOAs) of the pulsar p due to its generally unseen orbiting com- 
panion c. Its Keplerian expression is |27tl28] 



(cos E — e) sin oj + yl — e 2 sin E cos u 



(24) 



where E is the eccentric anomaly, and x p = a p sinl/c is the projected semi- 
major axis of the the pulsar's barycentric orbit; if m p = M, m c = m, it is 
Bp = om = mM^a. 

In general, the variation Ay of an observable Y over a full orbital revo- 
lution can be computed as 



AY 



dY . , 



2n 



dY dE dM ^ dY dtp 



dEdM dt 



dt , , 
lE UE - 



(25) 

where Ai is the mean anomaly and tp collectively denotes the other Keple- 
rian orbital elements. The rates Ad, if) entering eq. (|25p are instantaneous, 
and are obtained by computing the right-hand-sides of, e.g., the Lagrange 
equations onto the unperturbed Keplerian ellipse without averaging them 
over Pb- 

By using eq. (I25p . it turns out that the post-Keplerian time shift per- 
turbation due to MOND is 

7aA Pg /m/i b \ ( a \ 3 r ^f^e 2 \ . 

Ar D = — 7T — ev 1 — e z 1H cos a; smi. (26) 

p 8vrc V^b / V^m; V 2^ K ' 



2.3 The radial velocity 

The radial velocity V p [29] is a standard observable in spectroscopic studies 
of binaries [30]. Up to the radial velocity of the binary's center of mass Vo, 
the Keplerian expression of the radial velocity of the visible component of 
the binary reads 

V p = K [e cos uj + cos (oj + /)] = b lc ' — [e cos oj + cos (lo + /)] , (27) 

yl-e 2 

where a\ c is the barycentric semimajor axis of the binary's component whose 
light curve is available; K is the semi-amplitude of the radial velocity. In the 
case of extrasolar planetary systems, it is usually the hosting star; thus, a\ c = 
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a,M = mM^ a. In the case of spectroscopic binary stars, it may happen that 
the light curves of both the components (double lined spectroscopic binary 
stars) are available. The MOND-induced post-Keplerian perturbation of eq. 
(|2T|) turns out to be 

AV P =^A^(^)(l)\(l + ^)s^. (28) 

3 Confrontation with the observations 
3.1 Planets of the Solar System 

As far as the Solar System is concerned, tjvi = 39 d; thus the quasi-staticity 
condition is fully satisfied by the gaseous giant planets for which it is Pb ^ 
4300 d. Among them, Saturn, whose period amounts to Pb = 10759 d, 
is the most suitable to effectively constrain a since its orbit is nowadays 
known with ~ 20 m accuracy |25j thanks to the multi-year record of ac- 
curate radio-technical data from the Cassini spacecraft. In terms of its 
perihelion, its anomalous precession can nowadays be constrained down to 
sub-milliarcseconds per century (mas cty -1 ) level, as shown by Table [TJ If 

Table 1: Supplementary precessions A/2, Am of the longitudes of the node 
and of the perihelion for some planets of the Solar System estimated by 
Fienga et al. |25| with the INPOPlOa ephemerides. Data from Messenger 
and Cassini were used. Fienga et al. [25] fully modeled all standard Newto- 
nian/Einsteinian dynamics, apart from the solar Lense-Thirring effect which, 
however, is relevant only for Mercury. The reference {x, y} plane is the mean 
Earth's equator at J2000.0. The units are milliarcseconds per century (mas 
cty" 1 ). 





AQ (mas cty 1 ) 


Aw (mas cty 1 ) 


Mercury 


1.4 ± 1.8 


0.4 ±0.6 


Venus 


0.2 ± 1.5 


0.2 ± 1.5 


Earth 


0.0 ±0.9 


-0.2 ±0.9 


Mars 


-0.05 ±0.13 


-0.04 ±0.15 


Saturn 


-0.1 ±0.4 


0.15 ±0.65 



the case a K = n 2 a K= i, with a K= % ~ 1 is considered, eq. ([23]) and Table [T] 



S 



yield 



k\ < 2.5 x 10 5 . 



(29) 



Our result supports the guess by Milgrom [T3] that values of k > 10 5 might 
be excluded. The outer planets (Uranus, Neptune, Pluto) are not yet suit- 
able for such kind of analyses: indeed, their orbits are still poorly known 
because of a lack of extended records of radio-technical data. As far as their 
perihelia are concerned, their anomalous precessions are constrained to a 
4 — 5 arcseconds per century (" cty -1 ) level |31| . To be more precise, from 
eq. (f23j) for moderate eccentricities it turns out 



In addition to Saturn (m = 5.7 x 10 26 kg, a = 9.5 au), let us consider Pluto 
(m = 1.3 x 10 22 kg, a = 39.2 au); Pitjeva [S] yields Aw = 2.84 ± 4.51 " 
cty -1 for its anomalous perihelion precession. Thus, 



the constraint on k from Pluto would be 1460 times less tight than eq. 
(|29p inferred from Saturn. Although the orbit determination of Pluto will 
be improved by the ongoing New Horizons mission [32j to its system, its 
perihelion precession should be constrained down to a totally unrealistic 
0.001 mas cty -1 level in order to yield constraints competitive with eq. 
(f29j) . An analogous calculation for Neptune (m = lx 10 26 kg, a = 30.1 au, 
Am = —4.44 ±5.40 " cty -1 [31]) yields KNeptune ~ 28«;s a tum- It implies that 
the anomalous perihelion precession of Neptune should be improved down 
to a 0.1 mas cty -1 level. At present, no missions to the Neptunian system 
are scheduled. Nonetheless, the OSS (Outer Solar System) mission [33], 
aimed to test fundamental and planetary physics with Neptune, Triton and 
the Kuiper Belt, has been recently proposed; further studies are required to 
investigate the possibility that, as a potential by-product of OSS, the orbit 
determination of Neptune can reach the aforementioned demanding level 
of accuracy. The situation for Jupiter (m = 1.898 x 10 27 kg, a = 5.2 au) 
is, in perspective, more promising. At present, its perihelion precession is 
modestly constrained at a —41 ± 42 mas cty -1 level |25j : thus it is currently 
not competitive with Saturn. A 0.1 mas cty -1 level would be required: such 
a goal may, perhaps, not be too unrealistic in view of the ongoing Juno 
mission [34] . which should reach Jupiter in 2016 for a year-long scientific 




(30) 




(31) 
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phase, and of the approved_| JUICE mission [35] > to be launched in 2022, 
whose expected lifetime in the Jovian system is 3.5 yr. 

3.2 Radial velocities in binaries 

In general, according to eq. (|28p . the most potentially promising binaries 
are necessarily those orbiting slowly enough to fulfil the quasi-staticity con- 
dition. Moreover, they should move in highly elliptical, non-face-on orbits, 
and their masses should be comparable. Finally, the data records should 
cover at least one full orbital revolution. 

3.2.1 Exoplanets 

The wealth of exoplanets discovered so far allows, at least in principle, to 
select some of them for our purposes. 

Let us consider 55 Cnc d [36] which is a Jupiter-like planet (msin/ = 
3.835mj) orbiting a Sun-like star (M = 0.94M Q ; £m = 38 d) along a mod- 
erately elliptic orbit (e = 0.025) with a period P h = 14.28 yr = 5218 d; 
the other relevant parameters are u = 181.3°,/ = 53°. It was discovered 
spectroscopically; the accuracy in measuring the amplitude K of its radial 
velocity is [36] 

a K = 1.8 m s -1 . (32) 
By using eq. (|28j) for 55 Cnc d and eq. ([32]) . it turns out 

M<7xl0 8 , (33) 

which is 3 orders of magnitude weaker than the constraint of eq. (J29]) inferred 
from the perihelion precession of Saturn. It should be noticed that the use 
of eq. (I28p . which refers to the shift of the radial velocity over one full 
orbital revolution, is fully justified since Fischer et al. [36J analyzed 18 years 
of Doppler shift measurements of 55 Cnc. 

Other wide systems may yield better constraints, although not yet com- 
petitive with those from our Solar System. For example, HD 168443c [37] 
(M = O.995M , msin/ = 17.193mj, t M = 39.8 d, P b = 4.79 yr = 1749.83 
d, e = 0.2113, a = 2.8373 au, u = 64.87°, a K = 0.68 m s" 1 ) yields 

\k\ < 3 x 10 7 (34) 

by assuming / = 50° . Also in this case the use of eq. (f28l) is justified since the 
spectroscopic Doppler measurements cover more than one orbital period. A 

E http : //sci . esa. int/ science- e/www/object /index. cfm?f object id=50321 
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similar result may occur for 47 Uma d [38] (M = 1.03Mq, msinl = 1.6mj, 
t M = 40 d, P b = 38.3 yr = 14002 d, e = 0.16, a = 11.6 au, u = 110°, 
ok = 2.9 m s _1 ), but, in this case, the data used by Gregory et al. [38] span 
a period of just 21.6 years. 

3.2.2 Spectroscopic stellar binaries 

Looking at double lined spectroscopic binary stars, an interesting candidate 
is the a Cen AB system [39]. It is constituted by two Sun-like main se- 
quence stars A (M = I.IOSMq) and B (m = O.934M ) revolving along a 
wide (a = 23.52 au) and eccentric (e = 0.5179) orbit with Pb = 79.91 yr = 
29187.12 d » tM = 56.35 d. The standard deviations of their radial veloci- 
ties are [39] a ja) = 7.6 m s -1 , a (b) = 4.3 m s . Thus, from eq. (|28j) we 
obtain the tight constraints 



Such bounds are one order of magnitude tighter than the one of eq. (I29p 
inferred from the perihelion precession of Saturn. 

Other aspects of MOND, different from the effect treated here, were 
investigated with Proxima Centaur]! (a Cen C) [41H43] . 

3.3 Pulsars 

In order to fruitfully use eq. (1260 . the orbital period of the binary chosen 
should be larger than tM 46.7 d, obtained by using the standard value for 
the pulsar's mass M = 1.4M ; this implies that wide orbits are required. 
Moreover, they should be rather eccentric as well, and the mass m of the 
companion should not be too small with respect to the pulsar's one. Finally, 
timing observations should cover at least one full orbital revolution. As 
a consequence, most of the currently known binaries hosting at least one 
radiopulsar are to be excluded because they are often tight systems with 
very short periods. 



A partial exception is represented by the Earth-like planets |44j C (Pb = 
66.5 d, m = 0.0163mj, a = 0.36 au, e = 0.0186, / = 53°, u = 250.4) and D 
(P b = 98.2 d, m = 0.0164mj, a = 0.46 au, e = 0.0252, / = 47°, u = 108.3) 

6 It should be gravitationally associated with a Cen AB [40]. Their mutual separation 
is 15,000 au [2D]. 



k\ < 6.2 x 10 4 (A), 



(35) 



k\ < 4.2 x 10 4 (B). 



(36) 



11 



discovered in 1991 around the PSR 1257+12 pulsar (M = 1.4M ) [35]; the 
post-fit residuals for the TOAs was a Tp = 3.0 /xs [SJ. Applying eq. ([26j) to 
D yields 



Such a constraints is far not competitive with those inferred from Saturn 
(Section El]) and a Cen AB (Section [3X5} . 

4 Summary and conclusions 

We looked at the newly predicted quadrupolar MOND effect occurring in 



non-spherical, isolated and quasi-static yP^ 3> tu = y GMtot^o c 2 J sys- 
tems in deep Newtonian regime, and calculated some orbital effects for a 
localized binary system in the framework of the QUMOND theory. 

In particular, we worked out the secular precession of the pericenter, and 
the radial velocity and timing shifts per revolution. Our results are exact in 
the sense that no simplifying assumptions about the orbital geometry were 
used. 

They all contain a multiplicative scaling parameter a which, in certain 
MOND theories encompassing another characteristic constant in addition 
to the usual MOND acceleration scale Aq, can be expressed as a K = K 2 a±, 
with k > l,cti ~ 1. 

By using the latest orbital determinations of the planet of the Solar Sys- 
tem, we inferred |k| < 2.5 x 10 5 from the supplementary precession of the 
perihelion of Saturn. The double lined spectroscopic binary a Cen AB al- 
lowed to obtain | n\ < 6.2 x 10 4 (A) , | n\ < 4.2 x 10 4 (B) from our prediction for 
the shift in the radial velocity. The bounds that can be obtained by extraso- 
lar planets, including also those orbiting pulsars, are not yet competitive. In 
general, the best candidates are binary systems made of comparable masses 
moving along accurately determined wide and highly eccentric orbits. 

Our constraints are to be intended as somewhat preliminary because, 
strictly speaking, they did not come from a targeted analysis in which the 
MOND dynamics was explicitly modeled in processing the observations and 
a dedicated solve-for MOND parameter such as k was determined along with 
the other ones. Nonetheless, they are useful as indicative of the potential- 
ity offered by the systems considered, and may focus the attention just to 
them for more refined analyses. Moreover, in interpreting our constraints 
on k, it must be recalled that they were obtained within the framework of 



k\ < 1.5 x 10 



12 



(37) 
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QUMOND, which is not the nonrelativistic limit of theories such as TeVeS 
predicting high values of k. 
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